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Abstract 

We introduce models of heterogeneous systems with finite connectivity defined on ran- 
dom graphs to capture finite-coordination effects on the low-temperature behavior of finite 
dimensional systems. Our models use a description in terms of small deviations of particle 
coordinates from a set of reference positions, particularly appropriate for the description of 
low-temperature phenomena. A Born-von-Karman type expansion with random coefficients 
is used to model effects of frozen heterogeneities. The key quantity appearing in the theo- 
retical description is a full distribution of effective single-site potentials which needs to be 
determined self-consistently. If microscopic interactions are harmonic, the effective single-site 
potentials turn out to be harmonic as well, and the distribution of these single-site poten- 
tials is equivalent to a distribution of localization lengths used earlier in the description of 
chemical gels. For structural glasses characterized by frustration and anharmonicities in the 
microscopic interactions, the distribution of single-site potentials involves anharmonicities 
of all orders, and both single-well and double well potentials are observed, the latter with 
a broad spectrum of barrier heights. The appearance of glassy phases at low temperatures 
is marked by the appearance of asymmetries in the distribution of single-site potentials, as 
previously observed for fully connected systems. Double-well potentials with a broad spec- 
trum of barrier heights and asymmetries would give rise to the well known universal glassy 
low temperature anomalies when quantum effects are taken into account. 
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1 Introduction 



Understanding the behaviour of amorphous systems, both at elevated and and at low tempera- 
tures, continues to be a challenging problem. 

Structural glasses, for instance, show intriguing dynamic and thermodynamic behaviour when 
they are formed, e.g., by cooling from the liquid state [1]. A very sharp rise of relaxation times 
from microscopic to macroscopically large within a fairly narrow temperature interval entails 
that such systems will fall out of (local) equilibrium at a (cooling-rate dependent) calorimetric 
glass-transition temperature T g , accompanied by a sharp drop of the specific heat c p (T) and 
thermal expansion coefficient a(T). Signatures of a dynamic singularity are, however, observed 
in many glasses at a temperature T c larger than that of the calorimetric transition (UEH3J. Such a 
dynamic singularity was shown to emerge within a mode coupling theory (MCT) — a description 
of the dynamics of supercooled liquids in terms of correlation functions — via non-linear feedback 
mechanisms and memory-effects. The dynamic singularity occurs at a temperature where the 
viscosity is still fairly low; for a detailed review, see [2]. Conversely, if the behaviour of the entropy 
of the supercooled liquid is extrapolated to lower temperatures, it would become smaller than 
the entropy of the crystal at the Kauzmann temperature Tk, which is speculated to mark the 
location of a hypothetical true equilibrium glass transition, as it would be observed in the limit of 
infinitely slow cooling |4j. The Kauzmann temperature is usually quite close to the temperature 
To marking a hypothetical divergence of the viscosity when fitted according to a Vogel-Fulcher- 
Tammann form, and the two temperatures are often identified (and ought to be according to 
classical semi-microscopic [5J and phenomenological [6] approaches to the glass transition). 

In a series of seminal papers, Kirkpatrick, Thirumalai, and Wolynes [H [9] pointed out a 
connection between the transition to frozen order in certain spin-glass models and the structural 
glass transition. The spin-glass models in question exhibit an equilibrium phase transition to a 
low-temperature phase with one 1 step replica-symmetry broken structure of the matrix of spin- 
glass oder parameters. The transition is discontinuous in the order parameters, but exhibits no 
latent heat. Remarkably there is also a dynamical freezing transition at a temperature that is 
different (and higher) than the temperature To at which the equilibrium transition occurs 0] - 
The phenomenology is strikingly similar to that found in structural glasses, where the dynamical 
transition is associated with the mode coupling transition, while the equilibrium transition would 
be associated with a hypothetical equilibrium transition at the Kauzmann temperature T/~, and 
in both cases, the equilibrium transition appears to be driven by an entropy crisis. The analogy 
was further elucidated and strengthened by studies of spherical version of the p-spin interaction 
spin-glass model \10\ [TT] , for which macroscopic dynamical evolution equations for correlation 
and response functions can be formulated in closed form and have in fact the structure of 
idealized MCT equations [11] . In this case - unlike in structural glasses — the off-equilibrium 
dynamics in the non-ergodic low-temperature phase is also understood in rather great detail 
PI- 
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There is, of course, a sense in which the analogy between the physics of spin-glasses in the 
p-spin class and the physics of structural glasses is purely formal, as spin-glasses are systems 
with quenched disorder describing the interaction of degrees of freedom that are localized at all 
temperatures, whereas in structural glasses the disorder is self-induced and localization occurs 
only in the low-temperature phase. The issue of quenched disorder can be addressed at least 
schematically by noting that models with quenched disorder can, and have been used as bona- 
fide approximants to describe the physics of models with self-induced disorder |T3l EJ [15] . 
yet the issue of localization remains. Significant progress has, however, been made towards a 
first principles description of glassy thermodynamics by Mezard and Parisi [16J, based on an 
investigation of weakly coupled real replicas of liquid systems, elaborating on and extending 
related ideas of Monasson |17j . Franz and Parisi \18\ 119] . A more recent development are so- 
called lattice glasses [201 ED [221 [221 12H 125]' *- e - lattice gas models with discrete degrees of 
freedom on (generalized) Bethe lattices. These models are neither locally disordered nor of 
infinite connectivity, allowing in this way for an ordered, crystalline state, and showing again 
glassy behavior due to self-induced disorder. 

On a quite different note, structural glasses and other amorphous systems were found to exhibit 
thermodynamic, dynamic and transport properties at low temperatures, which are strikingly dif- 
ferent from those of their crystalline counterparts [26J. At temperatures below 1 K, the specific 
heat of glasses varies approximately linearly with temperature T, the behaviour of the thermal 
conductivity is close to quadratic in T, whereas in crystals both quantities are cubic in temper- 
ature. The behaviour of acoustic and dielectric response is also anomalous in this temperature 
range, when compared with the corresponding behaviour of crystals |27| . 

At temperatures at around 10 K there is an anomaly of a different kind in the specific heat which 
appears as a bump in the ratio C(T)/T 3 — the so-called Bose peak. The thermal conductivity 
exhibits a plateau in the Bose peak region, and begins to rise again as the temperature is further 
increased. 

Below 1 K, the low-temperature anomalies, as the phenomena just described have become to 
be referred to, are attributed to the existence of tunnelling systems which are absent in ideal 
crystals [281 129] . Tunnelling systems are thought to be realized as double- well configurations 
in the potential energy landscape, between which atoms, or groups of atoms can tunnel, even 
when thermal energies to overcome the intervening potential energy barriers are not available. 
Double-well potentials with a broad range of asymmetries and barrier heights can give rise to 
the broad spectrum of low-energy tunneling-excitations responsible for the anomalous thermal 
properties, whereas their resonant and relaxational interaction with the phonon bath creates 
anomalous response properties. At temperatures (energies) in the Bose peak region, the main 
excitations appear to be soft quasi-harmonic. Their origin and precise nature is much debated 
(e.g. [301 ED E21 [331 EI] ) ; in any case, strong hybridization between these states and the phonon 
bath is likely to be an important feature. 



2 



One of the remarkable aspects of glassy low-temperature anomalies is their considerable degree 
of universality. Two main contending theorys are available to explain this fact: Following ideas 
of Yu and Leggett [35], universality has been understood as a collective effect due to interac- 
tions between quantized excitations, i.e. the tunnelling systems [36]; alternatively universality is 
regarded as a consequence of the irregularities of the potential energy landscape that are frozen 
in at the glass transition, leading to a broad spectrum of quantized low-energy tunneling excita- 
tions [371 138j . The potential energy landscape being itself a collective affair generated through 
the interactions of constituent particles, one is led to expect a certain degree of insensitivity to 
detail of the low-temperature anomalies. This line of reasoning has been made more explicit in 
[39], and the key feature responsible for the universality of the low-temperature anomalies was 
demonstrated to be a clear separation of the energy scales frozen in at the glass transition on the 
one hand side, and that responsible for the low temperature physics on the other hand side. A 
simple explanation for the mysterious so-called quantitative universality, according to which not 
only power-laws describing the variation of some thermodynamic functions with temperature, 
but also prefactors in these power-laws are insensitive to detail, was also provided [39] . 

On the other hand, the description of the glassy freezing transition in these models is not as one 
would want to have it for structural glasses, as the transition to the non-ergodic low-temperature 
phase is not of the discontinuous 1-step replica-symmetry breaking (RSB) type as in the p-spin 
models, so the description chosen in [371 133 [39] does appear to miss essential aspects of glass 
transition physics. To be fair though, it should be mentioned that those models were initially 
constructed mainly to provide a description of low-temperature phenomena. However, being 
microscopic the question as to their high-temperature properties naturally arises, and in this 
regime one might regard them as inadequate in their present form. Conversely, regarding low- 
temperature phenomena, it is difficult to conceive how models of the p-spin type would describe 
the physics of tunnelling systems; in their spherical mean-field variant in particular [lOUllj . they 
would give rise to harmonic effective single-site potentials which could not support tunnelling 
systems. It is thus perhaps fair to say that a description of both, the low-temperature phenomena, 
and glass-transition physics within a single unified approach still eludes us. 

The main purpose of the present investigation to improve upon the models investigated in 
[HTJ [381 EH] , an d consider variants in which interactions are finitely coordinated, keeping them 
random though, as a schematic way of modelling a frozen glassy state. Our investigation will 
exploit the fact that that the understanding of dilute random systems — introduced many years 
ago by Viana and Bray [3D] — has witnessed significant progress in the last few years, regarding 
their equilibrium properties [HI E21 E3] , as wen as their dynamic behaviour [31 1 1^ 1 B7 [ 
This progress was also the basis for the previously discussed lattice-glass models, which, however, 
due to their discrete nature do not allow for the description of low-temperature excitations. 

It is hoped that the additional element of realism introduced by considering finitely coordinated 
systems, respecting in this feature an important property of finite dimensional systems, could 
help bring those models closer to reality. 
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The remainder of the paper is organized as follows. In Sect. [2] we introduce our model, and 
demonstrate how a finitely coordinated random model arises as a natural candidate to describe 
the physics of a glassy low temperature phase. We will look in detail only at a simplified version 
with scalar degrees of freedom. We describe the solution of the model using replica and mean field 
techniques in Sect. However, in Sect. 13.11 we evaluate only a replica-symmetric (RS) approxi- 
mation to the full theory. It turns out to have a similar complexity as a 1-RSB solution in models 
with discrete degrees of freedom. The key quantity appearing in the theoretical description is a 
full distribution of effective single-site potentials which needs to be determined self-consistently. 
A population-based stochastic algorithm is used to solve these RS self-consistency equations. 
The general set-up is very flexible and can be investigated for a broad spectrums of assumptions 
about the nature of the microscopic interactions as well as the connectivity distributions. 

In Sect. H]we evaluate the theory for models in which the microscopic interactions are purely 
harmonic, in which case the effective single-site potentials turn out to be harmonic as well; 
the distribution of the single-site potentials is then equivalent to a distribution of localization 
lengths used earlier in the description of systems of such type [321 EH ED E2J E3] , which may 
be thought of as systems in their gel-phase. For structural glasses characterized by frustration 
and anharmonicities in the microscopic interactions, investigated in Sect [5j the distribution of 
single-site potentials involves anharmonicities of all orders, and both single-well and double well 
potentials are observed, the latter with a broad spectrum of barrier heights. Although initially 
designed as models for the description of the low-temperature phase, it makes perfect sense to 
consider them also at high temperatures, where the microscopic symmetries of the interaction 
energy are unbroken. The appearance of glassy phases, as the temperature is lowered is marked 
by the appearance of asymmetries in the distribution of single-site potentials, as previously 
observed for fully connected systems. We investigate the phase diagram for a model system 
of this kind and compute thermodynamic function across the transition. As in fully connected 
systems, the transition is found to be continuous in our anharmonic model, and in this aspect the 
model we are considering here remains deficient. We discuss possible remedies of this deficiency 
in a concluding section. 

Details of some of the calculations were relegated to appendices. The paper closes with a sum- 
mary and an outlook on future interesting lines of research. 



2 The Model 



We consider a many-particle system described by the Hamiltonian: 





4 



with an interaction energy taking the form 



(2) 



Our initial goal is to attempt a microscopic, albeit schematic description of the physics of an 
amorphous random solid. We proceed as follows. To begin with, we decompose each particle 
coordinate ri into a fixed reference coordinate Ri and a relative coordinate U{ characterizing 
the (small) excursion from their respective reference position. Supposing that the initial pair 
potential has an effectively finite range, the interaction potential can then be formulated in 
terms of local connectivity variables and a (random) potential depending only on the relative 
coordinates. 

The connectivity is characterized by a matrix C with elements Cij, such that Cjj = 1 for close 
particles, for which |J2j — Rj\ < r c , with r c denoting the interaction range of the initial pair 
potential <j>, and Cjj = if |rj — rj\ — \Ri — Rj\ 3> r c and so <p(ri — rj) ~ 0. The interaction 
potential ([2]) may then be rewritten as 



Here the symbol (i, j) is used to denote pairs of sites; non-zero pair potentials <j>ij and single-site 
potentials Vj may vary from pair to pair and site to site, thus introducing residual disorder 
(arising, e.g., from interactions between different species, random nonzero Ri — Rj, etc., in the 
case of pair potentials, and the coupling of different species to external fields in case of the 
single-site potentials); we will specify 4> and V only later. One may think of the (f>ij(ui — Uj) 
as of a Born-von Karman expansion of the interaction potential about the reference positions, 
carried to all orders. 

In case one is interested mainly in low-temperature phenomena, one would expect that some 
low order expansion of the faj should be sufficient to capture the essential physics. 

As a crucial last ingredient, which renders the model solvable, we consider connectivity matrices 
C = {cij} which are not created as a consequence of an underlying geometrical arrangement of 
the interacting particles as described above, but rather we take them to be random. We assume 
that the local coordination number Li of site i is distributed according to a distribution Pq(L), 
with finite average connectivity C = ^2 L LPc(L) = jjJ2i-Li, and that the probability that 
two sites % ^ j with coordination numbers Li, Lj are connected (i.e. that Cij = Cji = 1) is 
proportional to LiLj (no correlations). The adjacency matrices C in the ensemble determined 
by Pc(L) then have the following distribution 




(3) 



(id) 



P C (C|{L,}) = - UpcM 4,,.Ip,..v 



(4) 
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where N is a normalization constant. It iturns out that the distribution ^ is independent of 
the conditioning {L{\ in the large N limit, for any typical realization of the set {Li} compatible 
with a given distribution 

Pc(L) = ^Y, 5 ^ ( 5 ) 

i 

of coordination numbers. A priori single bond probabilities 

C ( C\ 

Pc(Cij) = + ( 1 - J *cy,0 • (6) 

compatible with the average connectivity C are also included in (j3J) , as a matter of convenience 
rather than necessity, as constraints coming from the full distribution of coordination numbers 
are also encoded in Other, equivalent formulations are possible as noted in Appendix El 
Note that for Pq{L) = ^|-exp(— C) the connectivity matrix is that of an Erdos-Renyi random 
graph. 

The residual randomness in the pair and single-site potentials is described by normalized prob- 
ability densities P^ [</>] and Py [V] . 

In its present formulation, the model can be solved within replica mean-field theory. The nature 
of the random pair and single-site potentials as well as the distribution of connectivities can be 
left open for the time being. Specifying them in different ways one has access to describing a 
variety of different physical systems. 



3 Replica Analysis 

The generator of all the relevant physical quantities is the free energy, —N(3f{(3) = ln(Zjv), and 
using the standard replica trick (log Z) = lim ra ^o n l log {Z n }, we need to calculate the average 
of the n-fold replicated partition function over the bond-disorder, i.e. the C, and the fy*: 

(z%) c ,<t, = ( / n d < ex p I - ( 13 e - «?) + e if) ( ? ) 

\ i,a { \ (i,j),a i,a ) ) I c ^ 

Free energies are manifestly self-averaging w.r.t. any form of (finite-dimensional) site disorder, 
which is why taking a ^-average is not required here. It is useful to introduce the short-hand 
notation u = [u , . . . ,w n ) for a replicated vector, and the replica interaction (bond) and single- 
site weights: 



U b (u,v,4>) = exp ^-pY^^-v^ , U s (u,V) = ex V (^-pY V ( ua )) 



(8) 
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The line of reasoning involved in performing the average over connectivity matrices for non- 
Poissonian connectivity distribution closely follows [53]. Constraints on the local connectivities 
in the C-average are enforced via the identity 

It is found that the average of the replicated partition function can be written as a functional 
integral over the 'replica density' 

i 

Enforcing its definition in terms of conjugate variables p{u) we obtain after some standard steps 
(for details see appendix A): 



(Zn)c,* = JfS VpVp 6XP { N (l( Gb[p] Gm[pJ] + Gs[p] ) } ' (U 



where 



G b [p] = l^j du dv p(u) Ub(u,v,(f)) p(v) 
G m [p,p } = J du dv p(u) p(u) (12) 
G s [p] = J2Pc(L) (hxJdup L {u)U s {u,V) 



The path integral is dominated by its saddle point, and functional variation with respect to 
p(u), and p(u) leads to the stationarity conditions: 

p{u) = cl^J dvU b (u,v,ct>) p(v)\ (13) 

(~\ V- n mr/ P L -\u)U 3 (u,V) \ 

P( u ) = Z^ Pc ^ L > L \ TT~ -L(-\ u F \r\ I ^ 14 
i \J dv P {v) U s {v,V)/ v 



3.1 Replica Symmetry 

In order to be able to take the limit n — > in (jlip -(|14p as required for the replica method, we 
need to make an ansatz for the replicated density p in (fTUj) and its conjugate p, which assumes 
certain symmetry properties under permutation of the replicas. Within the present paper we 
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shall only explicitly deal with a replica symmetric (RS) ansatz. Variants which would break the 
symmetry between replicas are fairly straightforward to write down (see Appendix but for 
the model class studied in the present paper they turn out to be extremely complex and difficult 
to handle numerically. The reader will appreciate this statement once the RS theory has been 
developed. 

The RS ansatz describes a situation with unbroken replica symmetry, and without loss of gener- 
ality both p(u) and its conjugate p(u) can be written as (functional) superpositions of products 
of single replica functions of the form 



i.e. the single replica functions in (|15p and (|16p are taken to have Gibbsian form, with functions 
ip and t/j denoting single-replica potentials that generate the Gibbs distributions in question. 
Note that Dtp and T>i/j are suitable measures in function space. With the prefactor C as chosen 
in (116j) . it turns out that, both ir and tt are normalized probability density functions (pdfs) 
over function space (see Appendix A). In terms of these specifications, the replicated partition 
function (jlip may be re-expressed as a functional integral over tt and tt in a way that allows to 
isolate the leading order in n (i.e. O(n) in the limit n — ► 0). One obtains an expression of the 
form 




(15) 



(16) 



Here we use the convention 




(17) 



(Zn)c^ ~ ThrVvexp[nN 




L 



}) , (18) 



On introducing 




dit dv exp 



(-/?(/(«) +g(v) + h(u,v))) 



(19) 



and short hands of the form Virlip] = Dtp Tr[ip] and {T)7t}l = Y[? f tt[^], one can express the 
three functionals appearing in (|18p as 




(20) 
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<?mb>", 71"] 



£>tt[V>] Vn[i>] In 



Z 



g s>L [Tf] ~ /{^ttIl (in 




(21) 



(22) 



All L (coordination number) summations are over non-negative integers (and the convention is 
used that an empty product evaluates to unity). 

Note that one can easily identify a bond (i.e. pair- interaction) term Qf,, a mixture term Q m , and 
coordination number L dependent single site terms Q s ,l- 

The precise form of the various functions and parameters is obtained from the stationarity 
condition with respect to variations of ir, which involves solving stationarity conditions w.r.t. 
the conjugate functional and 7} as well (Lagrange multipliers A, A are introduced to take care 
of the normalization constraint on ir and tt in the variational procedure). We finally obtain (for 
details see e.g. [HIES]) the following coupled set of integral equations for n and ir: 



TT 



TT 



L 
C 



(23) 
(24) 



in which the functions V] and i&[ip, <f>] are defined as 

L 



with 



Z^,4>{v) = j du exp ( - f3%l)(u) - (3<f>(u - vfj . 



(25) 
(26) 

(27) 



To summarize, we see that ir^ip] can be obtained from sampling from P^[(p} and ir[ift], while 
can be obtained from sampling from jjPc(L), Py[V], and This provides a strategy of 
solving the fixed point equations via a population dynamics based algorithm. 

In the saddle point the expression for the free energy (see Appendix [B]) simplifies to 



-Y, P c(L) J{Vtt} l (lnZ[*[{# z ,F]] 



(28) 
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The internal energy is given by: 



E{0) 



s.L 



where (-)b 2 denotes the Gibbs average corresponding to the bond energy 



(29) 



(30) 



while (-) s L is the Gibbs average corresponding to the single-site potential V]. Note that 

contributions to the internal energy arising from temperature-dependences of the distributions 
7r and tt vanish owing to the stationarity condition on /(/?) w.r.t ir and 7T. 

Thermal averages of arbitrary single-site observables are evaluated as Gibbs averages in the 
ensemble ([25]) of effective single-site potentials V e s = + which are distributed according 

to 

(31) 



P[Vos] = ^2 Pc{L) {Vtt} l U V eS -^f[{^} L ,V 



i.e. for an arbitrary observable of the form A = -k Y2i J ^-( u i) we have that its thermal average is 
given by 



VP[V C 



off J 



/ du A{u) exp(- (3V eS {u)) 



z[K 



cff J 



Note that for Poissonian random graphs, the ensemble of single-replica potentials ip is equivalent 
to the ensemble of effective single-site potentials V e R. This follows by comparison of (|24p and 
(I3ip on noting that ^Pc(L) = Pc(L — 1) for Poisson distributions. This implies in particular 
that the equivalence between the ensembles of single-replica potentials and effective single-site 
potentials is lost for systems with connectivity distributions other than Poissonian. 

We note that an alternative way of obtaining the replica symmetric theory is via the Bethe-Peierls 
iterative method for obtaining free energies on (locally) tree-like structures [32] . Appendix [C] 
gives a sketch of the reasoning for the present model class. 

This concludes the general framework. 
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3.2 Orthogonal Function Representation 

It should be clear that instead of working with the representation-free full functional set-up 
described above, we could have chosen to represent the single replica potentials ip and tp in 
terms of their expansions using a suitable complete set of basis functions {^j: 

^{u) =^2x IM ip l ,(u) = x- ip(u), (33) 

and analogously ^(u) = x^-(p{u). Integrals over the function space are then replaced by multiple 
integrals over expansion coefficients. It is convenient to choose the basis functions {^} to be 
orthonormal with respect to a scalar product, 

(<Pn I <Pv) = <W • ( 34 ) 
such that any function / spanned by the {9?^} has an expansion 

/(u) = /.£(u), /„ = (/| ¥7,). (35) 

Instead of (|15|) and (|16p one would write 

f exp\-f3x ■ (p(u a )} 

p(u) = J dir(x) J] (36) 

f t-t exp[— 8x ■ (p(u a )} 

ph = cjdm) n ( 37 ) 

for the replica density and its conjugate, with Z{x) and Z(x) fixed by a normalization condition. 
Shorthands such as dir(x) = dx tt(x) are being used in analogy to before. 

Repeating the above line of reasoning one would obtain 



tt(x) = J dTr(x) (6 (x-Mx^)))^ ( 38 ) 
<x) = Y,c Pc(L) J l[dn(x e ) (5[x-X L _ 1 ({x e }))) v (39) 



with 



L-l 



X L -i({2Li}) = bp\V) + ^2x i , (40) 

e 

X(x,<t>) = -r X {<p\ ^Z^) , (41) 
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and 

Zx,<t>{v) = j du exp (-fix ■ <p(u) - 0(f)(u - «)) (42) 

instead of f|23|) - (|2T|) . The distribution (|3ip of effective single site potentials translates in an 
obvious way to a corresponding distribution of their expansion coefficients in an orthogonal 
function representation, given by 

r L 

P(x cS ) = P c(L) / II^(^) ( S {^S-X L (m))) v • (43) 

L J I 

An obvious advantage of explicit representations of this form is that one may use it to formulate 
natural approximation schemes by truncating the expansions at some suitable finite order. A 
corresponding disadvantage would be that such truncations can be expected to be efficient only, 
if the system of basis functions is well adapted to the problem being studied. 



4 Harmonic Couplings 

The first model that we investigate is a random graph of harmonically coupled scalar degrees 
of freedom with a distribution of the coupling strengths. We assume that there is no single 
site potential, so that the system is translationally invariant. The interaction energy between 
connected vertices (i, j) of the graph is given by 

k- ■ 

<f>ij{ui -uj) = -?r(ui - uj) 2 (44) 

Models of this type have been used to describe gels; typically, uniform distributions of cross-link 
strengths P(kij) = 8(kij — k) where chosen in that context [491 [50l I51j . Alternatively models 
with random harmonic couplings on regular lattices [33] have been looked at in connection with 
Bose-peak phenomena. Note that [51] , which — unlike our present modelling — provides a 
description preserving information about spatial structure, has recently been used to provide a 
microscopic underpinning of a phenomenological model of an elasticity theory for soft random 
solids based on randomly varying elastic constants [53]. In the present investigation, we combine 
elements of random structure with random values for existing couplings. 

It turns out that the case of harmonically coupled degrees of freedom provides the only example 
where the population dynamics described above is closed within a parameter space of finite 
dimensionality. In all other cases, finite dimensional parameter spaces can only give variational 
approximations to a full solution. 

It should be noted, however, that harmonically coupled systems have unique ground states and 
therefore miss essential elements of glassy low-temperature physics. E.g. they are unable to 
support two-level tunneling systems in their potential energy landscape. 
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On the positive side, one would expect the replica symmetric approximation to be exact for such 
systems. 

Since the local variables can take any real value, we take the modified Hermite polynomials 
<Pim{ u ) = } Hg(u), with hn = \pK c l il \x\ as the appropriate set of basis functions, in which the 



H^(u) are Hermite polynomials satisfying the recursion relation 

H IM+1 (u) = 2uH IM (u)-2n H^u) , » = 0, 1, 2.., oo, H {u) = 1 . (45) 
The polynomials (p^ are orthonormal with respect to the scalar product: 

(f\g)= fdu exp(-u 2 ) f(u) g(u) (46) 



In the present case the family of orthogonal polynomials can be truncated after (f2, to leave the 
following set of basis functions: 

feW}= bs'vsr^sr ■ (47) 



In the fixed point equations (|39Mip . the integral in (f42l) . for <fi(u — v) = |(n — w) 2 , is Gaussian 
and can be done analytically, such that the are determined as: 



where 

-lo. 8X2 
a = 1 + 



.^2 

thus leaving us with algebraic updates for the population dynamics only. 

We observe that a positive X2 generates a positive X2, thus guaranteeing stability starting 
from distributions ir(x) and tc(x) which are non-zero for X2,^2 G 1R + only. We also note that 
the evolution of the xq and xq distribution does not couple back into the distribution of non- 
constant contribution to the ifj(u) and ip(u) respectively. Finally, it is interesting to observe that 
the distribution of xi, X2, X\ and X2, are temperature independent for the harmonic system, 
while the xo and xq distributions exhibit a non-trivial temperature dependence. 

The results are presented in the following figures. Fig [T] shows the distribution of the expansion 
coefficient X2 of the single-replica potential in an expansion in modified Hermite polynomials as 
explained above, for a Poissonian random graph. As mentioned before, this is in the present case 
equivalent to the distribution of the expansion coefficient x e ff,2 of the (purely harmonic) effective 



13 



single-site potential. The expansion coefficient x e ff,2 is in turn related to the localization length 
£ of the particles via 

r 2 = ^ . (49) 

The multiple peak structure in the distribution of localization lengths is notable, and related to 
the fact that localization lengths of different particles are to first approximation determined by 
their coordination — the clearly discernible 7 peaks in the left panel of Fig Q] corresponding to 
coordinations L = 1 (dangling bond) to L = 7. These results have in fact been obtained earlier 
via different routes [51]. A zoom into the low xi (large £) region reveals significant substructure 
which shows that correlation lengths of particles depend in fact also on coordinations of their 
neighbours (and on coordinations of next neighbours and so on). 




-1 1 2 3 4 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Figure 1: Distribution tt{x2) for a harmonically coupled system with homogeneous coupling 
strengths fey = 1 on a Poissonian random graph of average connectivity C = 4; the inverse 
temperature is (5 = 1. Left: full distribution; right: zoom into the low X2 region. The 5-peak at 
X2 = in the left panel is not drawn to scale; its total weight corresponds to the fraction of 
vertices not in the percolating cluster. 

In Fig[2]we show analogous results for a graph with non-Poissonian degree distribution, assuming 
connectivities of the form L = 2 + L' with V a Poissonian of mean 2. In effect, therefore, we 
do again have mean connectivity C = 4. However, there are no longer isolated clusters of finite 
size. 

For the non-Poissonian random graph, a comparison with the distribution of the expansion 
coefficient x e fr,2 of the effective single site potentials in Fig [3] reveals that these ensembles - 
while qualitatively similar — are indeed not identical. 
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As a last result for the harmonic system, we look at the distribution the expansion coefficient 
x eS ,2 f° r a system in which we have non-uniform coupling strengths in addition to locally varying 
connectivities. In the present case we took k%j uniformly distributed in [0,1], and choose a non- 
Poissonian graph structure of with C = 4 of the same type as before. Due to the variability in the 
kij, the peak structure in the x e ff ^-distribution is smeared, and it is no longer possible to associate 
a value of the localization length with a local connectivity, simply by looking at the distribution. 
However, the method of determining the distribution of effective single site potentials easily 
allows to monitor the contributions coming from sites with different connectivities, as illustrated 
in Fig. HI 



5 Glassy Sytems with Anharmonicities 

The next model that we investigate is a random graph of coupled (scalar) degrees of free- 
dom with a distribution P(kij) of harmonic coupling strengths and fixed anharmonic coupling 
strength. Furthermore, we assume that there are non-random anharmonic single-site potentials. 
The interaction energy between connected vertices of the graph is taken to be of the form 

k- ■ 

(j>ij{ui - uj) = -^-{v,i - Ujf + \ 4 {ui - Uj) 4 (50) 




Figure 2: Distribution tt{x2) for a harmonically coupled system with homogeneous coupling 
strengths k^ = 1 on a non-Poissonian random graph of average connectivity C = 4 as described 
in the text; the inverse temperature is (3 = 1. Left: full distribution; right: zoom into the low %i 
region. 
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x eff,2 

Figure 3: Distribution -P(x e ff,2) for a harmonically coupled system with homogeneous coupling 
strengths kij = 1 on a non-Poissonian random graph of average connectivity C = 4 as described 
in the text; the inverse temperature is (3 = 1. 
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x eff,2 

Figure 4: Distribution ^(2^,2) shown together with separate contributions to it coming from 
sites with connectivities L = 2,3, ... ,7, for a harmonically coupled system with couplings kij ~ 
U[0, 1] on a non-Poissonian random graph of average connectivity C = 4 as described in the 
text; the inverse temperature is (3 = 1. 

and we introduce a vertex-independent single-site potential 

V{ Ui ) = g 4 uf (51) 

If the support of P(kij) contains negative couplings, at least one of the quartic couplings A4 or 
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54 must be positive for the system to be stable. 




A fully connected variant of this model with 54 > and A4 = was studied in [38] , the 

corresponding model with 54 = and A4 > in For the purposes of describing glassy 
low-temperature anomalies either form of stabilization by a quartic anharmonicity appeared to 
be acceptable. Differences between low-temperature properties of the two models were mainly in 
the details, although from a fundamental point of view the system with 54 = would of course 
be preferable, since 54 7^ breaks translational invariance. 

For the purposes of the present investigation we shall nevertheless stick to the local stabilizing 
potentials, i.e. choose A4 = and 54 7^ 0, because this form of stabilization is numerically much 
easier to handle than its translationally invariant counterpart. We will turn to anharmonic 
systems with full translational invariance in a separate publication. 

From the solutions for the distributions 7r[V>] and fc[ip] of single-replica potentials one obtains a 
distribution of effective-single site potentials via (|31l) . 

Using an expansion of effective single-site potentials V e s{u) in terms of the modified Hermite 
polynomials used in Sec H] to describe harmonic systems, 



the distribution of effective-single site potentials translates into a distribution of expansion 
coefficients {x^} (to simplify notation, we drop the subscript 'eff' on the expansion coefficients). 

For systems with (quartic) anharmonicities, the solution of the fixed point equations (|23p . (|24p 
involves distributions of single-replica potentials i\) and ip with support on potentials having 
anharmonicities of all orders. Any finitely truncated orthogonal function representation of these 
distributions would therefore only generate an approximate solution to the full problem (on 
top of the approximation implied in the RS assumption). It is therefore advisable to use the full 
functional approach (not least in order to check the quality of finite dimensional approximations) . 
Below we therefore show distributions of some low order expansion coefficients as obtained in 
the functional approach and compare them with corresponding distributions obtained within a 
low-order truncation scheme. 

We fix the energy scale of the system by setting (74 = 1, and introduce frustration in the system 
by taking Gaussian distributed harmonic couplings, fcy ~ AA(0, <7fc). A system of this type then 
undergoes a transition from a symmetric high temperature phase, in which all effective single-site 
potentials respect the microscopic inversion symmetry of the original interaction energy V e s(u) = 
V e s (— u), to a glassy low-temperature phase, in which this symmetry is spontaneously broken, i.e. 
which is described by an ensemble of single-site potentials for which typically V c d(u) ^ V e g(—u). 
This transition occurs at a temperature T c which depends on the strength cr/t of the disorder 
(and of course on the connectivity distribution). Indeed, for a combination of quadratic and 




(52) 
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quartic potentials as chosen in the present setup, dimensional analysis reveals that the critical 
condition should have a parameter dependence of the form 

2 

T c (a k ,g 4 ) = ^T c (l,l) ■ (53) 

54 

as already found for the fully connected system [39]. Here the value of T c (l,l) still depends, 
of course, on the connectivity distribution. E.g., for the system with the non-Poissonian con- 
nectivity distribution with average degree C = 4 looked at earlier, we determined T c (l, 1) = 

0. 932 ± 0.005 /3 C (1, 1) — 1.073, and we verified the scaling (|53p with high precision within the 
full functional treatment of the self-consistency equations. It is important to note, however, that 
approximate treatments, such as the low-order orthogonal function approximation procedure, 
need not respect this scaling of T c (ak, 54), though we found them to be not far off. 

The transition to glassy order at low temperatures is continuous in the sense that widths of the 
distributions of the expansion coefficients x\, x$, X5, ... coupling to anti-symmetric functions go 
to zero continuously as the transition temperature is approached from below. Fig [5] illustrates 
this for a system with a non-Poissonian connectivity distribution of the type used in Sec 2] with 
C = 4, and = 1. If an orthogonal function representation of the single-replica functions 
in terms of Hermite polynomials, truncated at order four is used instead of the full functional 
approach, we find the order of the transition to be unchanged, though the transition point is 
shifted to a sligthly higher temperature, T c (l, 1) ~ 1.043 4=> /3 C (1, 1) — 0.959 (differing from the 
exact value by approximately 3%). 

We note that the population-dynamics algorithm exhibits critical slowing down as the transition 
is approached. In Fig[5j It is thererfore crucial to use sufficiently long runs and small temperature 
steps close to the transition to ensure that the algorithm has properly converged, before taking 
measurements. 

Fig [6] shows distributions of expansion coefficients for a system with the non-Poissonian connec- 
tivity distribution used before, at = 1.1, for T = 1 and T = 0, and compares distributions 
obtained from a full functional solution with those computed in an orthogonal function rep- 
resentation using Hermite polynomials that is truncated at order 4. Semilogarithmic plots are 
used to exhibit the behaviour in the tails of the distributions. At T = 1, it turns out that the 
distributions of odd coefficients xi, X3, . . . appear to be affected by the low order truncation in 
a slightly stronger way than the distributions of even oder coefficients, and distributions of low 
order odd coefficients appear to be systematically slightly wider when computed via a truncated 
expansion approach; the main reason seems to be that the low order truncation has slightly 
higher critical temperature and with identical parameter settings, the system appears to be 
deeper in the glassy phase when described within the low order truncation scheme than when 
described in terms of full functional self-consistency. The distributions of the even coefficients, 

1. e. X2, Xi ... computed in an approximation scheme truncated order 4 appear to be less strongly 
affected by the truncation scheme. 
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p 

Figure 5: Width a(x\) of the distribution of the first expansion coefficient x\ in in the ensemble 
of effective single-site potentials. The width goes to zero with a square-root singularity at (3 C ~ 
1.073. A fit displaying a square-root singularity with analytic corrections describes the data very 
well over a wide range of inverse temperatures. For this system = 1. 

It should be noted that there are strong correlations between the expansion coefficients that 
can be exhibited, e.g. by computing their correlation coefficients or by inspecting scatterplots 
of individual realizations. Fig. [7| shows some examples. For the system under investigation we 
found the (a^ifc-i) x 2l) scatterplots to be similar up to scaling and symmetry transformations; 
analogous similarities exist in (x2k, £2^)-f arm h es a t fixed k. 

Differences between the full functional solutions and those obtained within a low-order expansion 
become rapidly less pronounced when moving away from the transition temperature. This is 
particularly clear in the the zero-temperature limit where an orthogonal function representation 
truncated at order 4 describes the results remarkably well, as seen in the wider set of curves 
also displayed in Fig [6l This might have been anticipated, as anharmonicities typically play a 
less significant role at low temperatures. 

Another feature that deserves mention is that the distribution P{x2) of the second order coeffi- 
cient his hardly temperature dependent at all (e.g., the zero-temperature limit of this distribution 
is virtually indistinguishable from the one at T = 1 (as seen in the second panel of Fig [6]) whereas 
all other distributions do show a noticeable variation with temperature. 

Fig E] gives the phase diagram of the system, showing the transition temperature as a function 
of the strength of the disorder in the harmonic coupling constants. 

In Fig [9] we show free energy, internal energy, entropy and specific heat as functions of tem- 
perature, for <7fc = 1 and (74 = 1, so that T c ~ 0.93. The specific heat is obtained through a 
numerical differentiation from the internal energy, and becomes fairly noisy for temperatures 
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Figure 6: Distributions of low-order expansion coefficient x„, fi = 1, . . . ,8 in in the ensemble of 
effective single-site potentials for at = 1.1 at T = 1 (narrower set of curves) and T = (wider 
set of curves). Solutions obtained within the full functional approach (full lines) are compared 
with those obtained by truncating expansions of single-replica potentials in terms of Hermite 
polynomials at order 4. (dashed lines) of) 




Figure 7: Scatterplots showing 2 dimensional projections of vectors of expansion coefficients 
for realizations of effective potentials. First group of panels (xi, X2), (x±, X3), . . . (x±, x$), second 
group (x2, X4), (X2, x%), (X2, x%), third group (x4,X6) and (x±,x%), and last panel (xq,xs). 
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Figure 8: Transition temperature of the anharmonic system as a function of the standard de- 
viation cifc of the random harmonic coupling, and the strength of the quartic stabilizing 
potential. Iso-lines for T c = 0.1, T c = 1, T c = 10, T c = 100, and T c = 1000 are also shown in the 
(o-fc,54)-plane 

below T ~ 0.5. Notice that the entropy becomes negative for T < 0.31, i.e. well inside the glassy 
phase. As our system has continous degrees of freedom, no strong conclusions can be drawn from 
this observation, however. It must be pointed out, though, that we expect replica symmetry to 
be broken throughout the low-temperature phase, and that our results require further checks, 
e.g. through simulations. 

Finally, a word on numerics may perhaps be in order. It turns out that a good compromise 
between control and efficiency in evaluating functions like ^[■0,0] on a sufficiently dense grid 
of points v via ([27]) as required for the population dynamics algorithm, is to use Gaussian 
quadratures with sufficiently many grid points. Further acceleration of numerical procedures is 
possible by evaluating functions initially on coarser grids, and obtaining their values on the fine 
grids required for further processing via suitable interpolation algorithms. 
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Figure 9: Thermodynamic functions: free Energy / (first panel), internal energy E (second 
panel), entropy S (third panel), and specific heat C, (fourth panel). 

6 Conclusion 

We have introduced and solved a wide class of random models with finite-connectivity for the 
description of amorphous solid phases. The description is in terms of deviations of particle 
positions from a set of reference positions. The framework is kept general in the sense that we 
do not need not restrict the nature of the pair potential, nor the nature of the connectivity 
distribution, as long as they allow taking a meaningful thermodynamic limit. By varying the 
connectivity distributions or the nature of interaction potentials we expect to be able to describe, 
at least in a schematic manner, a fairly wide spectrum of different physical systems. 

In the present paper we have looked in detail only at harmonically coupled systems with scalar 
degrees of freedom, thought to provide an appropriate description of systems in their gel-phase, 
and at a model with local quartic anharmonicities which could be used to describe glassy low- 
temperature properties. The transition to the phase with frozen glassy order in the model con- 
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sidered here is continuous, as in the corresponding fully connected variant. In this respect, the 
description of glass-transition physics within this model must still be regarded as deficient. This 
may or may not be due to the fact that the low-order expansions of the interaction energy chosen 
in the present paper can be considered as adequate description of the physics only in the limit 
of low temperatures. 

One promising possibility to improve the description of glass-transition physics within the present 
model-class would be to look at models with non-confining interaction potentials. This would 
inter alia allow obtaining more adequate descriptions of high temperature phases, as phases 
with vanishing inverse localization lengths (delocalized particles). It remains to be seen whether 
this modification could put our model into the class for which glass-transition physics would 
be described by the succession of dynamic freezing and entropy-crisis scenario thought to be 
appropriate for structural glasses. 

Other problems we might be able to shed some analytic light on are correlations between co- 
ordination and effective potential statistics, or — by looking at mixtures of different particle 
species — the correlations between particle species and their enhanced or depressed role in the 
formation of tunneling centres in amorphous systems. 

Another open question is related to the fact that the present models — unlike their fully con- 
nected counterparts — are not completely described in terms of ensembles of effective non- 
interacting single-site problems. The two-particle contribution appearing in the free energy is a 
consequence of the fact that finite connectivity systems of the type considered here would no 
longer support simple single particle excitations. It is conceivable that the two-particle contri- 
bution would allow an interpretation similar to the 'heat-bath' of phonons acting on tunnelling 
systems that is introduced in an ad-hoc manner into phenomenological models of glassy low- 
temperature anomalies [28 | |29 | 155]. 

We intend to look into some of these questions in future publications. 

Acknowledgements It is a pleasure to thank ACC Coolen, J Hatchett, I Perez-Castillo, and 
P Sollich for many illuminating discussions and helpful remarks. 
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A Details of the Replica Analysis 



A.l Averaging the Replicated Partition Function 

Given the distribution @ of connectivity matrices, the average (|7J) of the replicated partition 
function has the following structure 



(7 n \ = — = - L2L (ka\ 

Note that both the Cij-summations and the average over local pair potentials (pij, factor w.r.t 
the pairs 

It should be mentioned that other, equivalent ways of introducing the ensemble of connectivity 
matrices are possible. E.g. one could include summations over all possible assignments of local 
connectivities {Li} compatible with a given Pc(L) in both T and J\f in (154p . without changing 
final results, as this would just create contributions coming from the entropy of the connectivity 
distribution in T and M, which cancel. Similarly, one could omit including the single bond 
distributions pc( c ij) m T an d M, as long as the constraints required by the desired connectivity 
contribution Pq(L) are kept — at the cost of considerably complicating the algebra, and without 
changing final results. In what follows, we stick to the version embodied in (|54H . 

To evaluate (|54p . one uses integral representations of the form ((9]) to express the constraints on 
the local coordinations Lj. The normalization constant N in (|54h is then transformed into an 
integral expression, which can be evaluated via the saddle point method, along the following 
transparent sequence of steps 



M = II Tr /'ri^M,..,,. II / 



d^i (Sj(^i) c ij— Li) 



2mz. 



11 ./ 2vriz ?: 1 



- /^-H^" 11 "^ ^? '"(«)}) (55) 

Here, we used large N asymptotics to re-exponentiate the Cjj averages in line 2, and introduced 
Po = ~k Yli z ii enforced by the conjugate po, to obtain an expression that can be evaluated by 
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the saddle point method, and finally exploited the fact that (|9|) entails 

(56) 



2ttiz L\ dz L 



2 = 



As the Li are distributed according to Pc(L), the last contribution in the exponent in (|55p can 
be expressed as the L-average ^2iPc{L) In (tj) by the law of large numbers. In the saddle 
point of (f55l) . i.e. (po = 1, /5o = C) we finally have: 

AT ~ exp fi\T P C (L) In ^ - C j j (57) 

Given any typical realization of local potentials {Vi,i = 1, .., N}, and using the shorthands © 
for the replicated bond and site weights we have for nominator T in (154j) . following similar lines 
of reasoning: 

T = FI / duiU s (ui,Vi) TTTr pc(c»i) ^c lJlCjl ((u b (ui,uj,<f)j) lJ ) Y[S Li ,Y: 

1 (m) 1 

y J 27TlZj 

z 

= y PpPp exp (iV - 1) - G m [p, Pi}) II/ ^) / 

= J VpVp exp ^N^(G b [p]-l)-G m [p,p} + G s [p]]^ (58) 

The definition (fT0|) of the replica density enforced by a conjugate p(tt) was used to express the 
result in terms of a path-integral that can be evaluated by the saddle point method. 

The three functionals appearing in (|58p are given as 

G b [p] = (^J du dv p(u) Ub(u,v,<t))p(v)j 
G m [p,p] = j du p(u) p(u) (59) 
G.\p] = ^P c (L)(lnjdu^-U s (u,V)^ (60) 

where the law of large numbers is used to write G s as an average over the site disorder (in the 
Li and the Vi). 



2mZi 1 
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A. 2 Replica Symmetry 



We proceed by making an ansatz for p and p which amounts to assuming unbroken replica 
symmetry (RS): 

exp[-/3V>(it a )] 



p{u) 
p(u) 



a 

t>4> m n 



z\ 



exp[-[3Tp(u a )\ 



(61) 
(62) 



where tt and n are functionals over the space of single site potentials and Z[tp] and Z[tp] are 
defined via (fT7l) . We then obtain: 

Z2[ipi,ip2,<f>] 



G b [p] ^ pl + n f Ptt^iIPtt^] An 
G m [p,p] ~ PoPo + n / Vtt^Vt!:^] In 



Z[fl>]Z$\) 



(63) 



G,[/5 



where the definition (I19p of Z2, the convention {"D7r}i = Y[f ^7r[^] , and the abbreviations 
Po = f Vir[ip] , and Po = f T>n[ip] for full integrals over the functionals 7r and n have been used. 

We first solve the saddle point equations to 0(1) = O(n ) in n, to obtain po = 1> and /?o = C. 
This entails that the O(n ) contributions to T cancel with J\f, and we finally obtain 



(Z n N ) ( 



J Vtt Pvrexp ^nJV j^&M - CO^vr] + ^P C (L) ^[tt] J \ 



in which 7r and tt are now normalised functionals, and where 



^ m [7r,7r] ~ / Vir[%[)} Vn[ip] In 
g s , L [vr] ~ l{Vn} L ( In ' Z 



Zty]Z[ 



(64) 

(65) 
(66) 

(67) 



The saddle point equations w.r.t. the normalised it, tt can then be expressed in the form ([23 
and (f24l) as given in Sect 13.11 
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A. 3 Breaking Replica Symmetry 



It is expected that replica symmetry will be broken in phases with frozen glassy order. We will in 
the present paper not go on to solve our models assuming phases with broken replica symmetry. 
However we would here like to present at least the ansatz which would describe phases with 
one level of replica symmetry breaking in the spirit of Parisi' hierarchical scheme. It is based 
on grouping the n replica in n/m blocks containing m replica each, and assuming symmetry of 
the solution w.r.t. permutations of replica within blocks though not between blocks. This would 
lead to 



exp[-^ fc (« a )] 



p(u) = fvir V[tt] f[ f V^ k Trtyfc] J] 

3 k=l J a=(fc-l)m+l ^ feJ 

r n / m r km r a 7 I aw 

m = \m vm n / v** m n eM ~ PMu )] 



a=(k— l)m+l 



(68) 
(69) 



as ansatz for the replica-density (llOh and the conjugate required to enforce its definition, re- 
spectively. I.e., we get a functional superposition of products of block-replica densities, each of 
which a functional superposition of products of single-replica Gibbs distributions over the replica 
within a block. 

This ansats would again allow to formulate the n — > 0-limit of the theory, and would lead to 
fixed point equations for the weight functions V, it, V and n that can be cast into a form which 
would admit a solution in terms of a stochastic population-based algorithm — an algorithm, 
however, which requires to maintain an ensemble of populations of the type used in the RS 
version of the theory. In addition, there is the stationarity requirement on the Free Energy with 
respect to the partitioning parameter m. The numerical effort required to solve this problem 
will be considerable. 



B Free Energy, Thermodynamic Functions 

The replica symmetric expression of the free energy as a funcional of tt and n obtained from 
(fT8 >(p22"]l initially gives 



C 



Vtt[iP]Vtt[iP'} (In 



C / £>7r[V>]£>7rk/'] In 



J2 P c(L) /{Pvt}l( In 



Z 



nil z 



(70) 
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in which n and tt are taken to be solutions of the fixed-point equations. 

Replacing logarithms of products (quotients) by sums (differences) of logaritms, and recalling 
and (L) = C, we first get 



C 



Pf(P) = -5- / T>ir[ip]T>ir[i//] (in [Z 2 [ip,ip f , 



C \ V7r[ip}VTr[i>] In Z[tp + ip] 



+ Y J P c{L) {Vtt} l (in 



(71) 



v 



In a second step, we express 7r[?/>] appearing in the second integral by an integral over Tr[ip] via 
(1231) to get 



£>7r[V>]£>7r[V>]ln Z[ifj + $ 



In 



V7c[ip]Vn[ip'] J Vip {5 ^ - V[ip', (f) 
Vir[ip}Vir[ip'} (in Z[ip + $[ip', <f> 
V^pTT^'} (In [Z 2 [^,0]]). 



In the last step, the definition (|26p . (|27p of ^[-i//, 0]) was invoked. We thus note that the first 
and second contribution to (|7ip can be combined so as to finally obtain (1281) . 

The internal energy is obtained from the free energy via the relation E(/3) = ^(J3f(J3)). 



C Bethe-Peierls Method 



The replica-symmetric fixed point equations can be derived as resulting from a recursive evalu- 
ation of the partition function of the system described by an energy function of the form 

= Cij<j>ij(ui - Uj) + V{ui) (72) 

assuming that the collection of nodes form a tree-like structure. Scalar degrees of freedom and 
site-independent on-site potentials V are assumed for simplicity. 

The restricted partition function computed at fixed value ui for the coordinate at i is 

Zi(ui) = ex.p{-(3V(ui)} Y{ duj e w{-P<l>ij('Ui-Uj)}Zj\i(uj) (73) 

in which denotes the set of vertices connected to vertex i, the Zju(uj) are the restricted 
partition functions of the subtrees rooted in j - exluding the tree rooted in node i - evaluated 
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at the value Uj of the coordinate at the root-node j. If the system were tree- like, the Uj integrals 
in this expression would factor, 

Zi(ui) = exp{-/3V(ui)} Y[ / dujexp{-/3(f>ij(ui - u^jZ^Uj), (74) 

and we would have the recursion 

Zj\i(uj) = exp{-pV(uj)} J| du k exp{-P4> jk (uj - u k )}Z klj (u k ) . (75) 

keM(j)\i 

This recursion is the basis for the iterative Bethe-Peierls method for solving statistical mechanical 
problems on tree-like structures. 

The system we are considering is, however, only locally tree-like, i.e. there are loops with lengths 
of O(lniV). Factorization is, therfore, in principle only approximate. However, if the system has 
only a single thermodynamic state and correlations decay exponentially with distance (this is 
the situation of unbroken replica symmetry), then factorization would be asymptotically exact 
in the thermodynamic limit, as common ancestors of two sites neighbouring on a given site 
occur typically at distance O(lniV) if the given site is removed, entailing that correlations are 
negligible in the thermodynamic limit. 

One introduces cavity potentials ^juiui) associated with the factors appearing in (|74"|) via 



exp{-/3Vj|i(iii)} = / dujexp{-(3cpij(ui - Uj)}Zj\i{uj) (76) 

For the restricted partition functions Zj^(uj) of the sub-trees themselves we have 

Z^iuj) = exp{-pV{uj)} Y\ ex P{-/%|«( u i)} ( 77 ) 

keAf(j)\i 

In terms of the cavity potentials ipju and the free energies ipju{v,j) associated with restricted 
sub-tree partition functions via 



z j\i( u j) = ex v{-/34>j\i(uj)} (78) 



one has the recursion 



^j\i( u j) = v ( u j) + ^2 (79) 

keAf(j)\i 

^k\j{uj) = - P' 1 In J du k exp{-(34> jk (uj - u k ) - Pip k ij(uj)} 

= -p-HnZ^^ (80) 
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The replica symmetric fixed-point equations derived in Sect. (|3.ip are now recovered by observing 
that due to the locally varying coordination and the heterogeneity of the interaction potentials 
4>ij associated with the edges (i,j) of the tree-like graph, the cavity potentials and the restricted 
sub-tree free energies are random. Hence the iteration has to be formulated in a probabilistic 
setting as an iteration for probabilities of finding certain restricted sub-tree free energies ip and 
cavity potentials if). For a homogeneously disordered (locally) tree-like graph this would lead to 
self-consistency conditions that have to hold in the thermodynamic limit. Using Eqs. (I79p and 
([BOD one finds 



l r L ~ l i 

*w = Ec Pc(L) /n p *H^" fi - i[{ ^ }1 l (81) 

L>0 J k=l 

m = fTht$\(8$ -%,</>]]) (82) 



for the self-consistency conditions, on a graph with connectivity distribution Pc(L) of average 
coordination C, and \P l-i and ^[^,(f) u ] are as defined via ()25[) and (|26p . 
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